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APPROXIMATE INDICIAL LIFT FUNCTIONS FOR SEVERAL WINGS 
OF FINITE SPAN IN INCOMPRESSIBLE FLOW AS OBTAINED 
FROM OSCILLATORY LEFT COEFFICIENTS 
By Joseph A. Drischler 

SUMMARY 


The unsteady-lift functions for a wing undergoing a sudden change 
in sinking speed have been presented for delta wings having aspect ratios 
of 0, 2, and !(- and for rectangular and elliptical wings having aspect 
ratios of 0, 3y and 6. For the elliptical and rectangular wings the 
spanwise lift distributions were also presented. The se functions were 
calculated from the lift coefficients associated with a wing oscillating 
harmonically in pure translational motion, as obtained from several 
sources. 

The results of these calculations indicate that the normalized 
unsteady-lift functions are substantially independent of the shape of 
the plan form for elliptical, rectangular, or moderately tapered wings; 
however, for delta wings the increase of lift toward the steady-state 
value is much more rapid than that for the aforementioned wings of the 
same aspect ratio. These results also corroborate the results of other 
investigations in that the rate of growth of lift tends to increase with 
a decrease in aspect ratio. The shape of the spanwise distributions of 
the indicial lift seems to be, for al 1 practical purposes, independent of 
time for rectangular and elliptical wings. 


INTRODUCTION 


The indicial lift functions due to a sudden change in sinking speed 
and due to the penetration of a sharp -edge normal gust (hereinafter des- 
ignated as k 2 _(s) and k 2 (s), respectively j are fundamental in the cal- 
culation of airplane transient motions and gust loads. These functions 
have been determined, however, for only a few wing configurations in 
incompres sible flow and for only wings of infinite aspect ratio in com- 
pressible subsonic flow. Initial derivations of the functions k-^Cs) 

and k 2 (s) were made for two-dimensional incompres sible flow by Wagner 
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(ref. 1) and Kussner (ref. 2), respectively; in reference 3 Von Ksrman 
and Sears pointed out and corrected an error in reference 2. Indicia! 
lift functions for incompressible flow have been calculated by Robert T. 
Jones (ref. 4) for elliptical wings having aspect ratios of 3 and 6 and 
by W. Prichard Jones (ref. 5) for rectangular wings with aspect ratios 
of 4 and 6 and also for a moderately tapered wing having an aspect ratio 
of For subsonic compressible flow, the indicial lift functions 

have been determined for infinite-aspect-ratio wings in references 6 , 7 , 
and 8 . No information appears to be available, however, on these func- 
tions for wings of finite aspect ratio in subsonic compressible flow, 
although current work on oscillatory lift coefficients may in time lead 
to further information on the indicial lift functions for this condition. 

The main purpose of the present report is to extend the knowledge 
of the incompressible -flow indicial function k^(s) to a plan form of 

current interest - namely, the delta wing - and to make a comparison of 
the k-j_ functions for elliptical, rectangular, tapered, and delta wings 

in order to assess the effects of plan-form shapes and aspect ratios. 

An additional purpose is to present some information on the variation 
of the spanwise distribution of the indicial lift function k^(s) since 

no published information on this distribution has been found. Such 
information is of interest in the determination of transient stresses in 
wing structures. 

The method used to obtain the k^ function was similar to that used 

to obtain the functions for a wing in two-dimensional subsonic compress- 
ible flow, as reported in references 6 and 7 * The indicial functions 
were obtained from their associated oscillatory lift coefficients (com- 
monly used in flutter and frequency- response calculations) by means of 
Garrick's reciprocal relations (ref. 9 ) • The oscillatory coefficients 
were determined by means of the methods presented by Lawrence and Gerber 
in reference 10 and by ReiBsner and Stevens in reference 11. (The method 
of ref. 11 also allows calculation of local oscillatory lift coefficients, 
leading to spanwise distribution of the indicial lift functions.) 

The indicial lift functions k^(s) were derived for delta wings 
with aspect ratios of 2 and b, and the spanwise variation of k^ functions 

was determined for elliptical and rectangular wings with aspect ratios 
of 3 and 6 . Also calculated were the unsteady- lift functions k-[_(s) 

and k 2 (s) for elliptical, rectangular, and delta wings of vanishingly 

small aspect ratio by a method based on an extension of the low-aspect- 
ratio theory of Robert T. Jones (ref. 12). The unsteady- lift functions 
previously obtained in references 1 to 5 are presented herein to show 
the effects of wing plan-form shape and aspect ratio. The limitations 
of the methods used in calculating the oscillatory coefficients precluded 
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the calculation of the spanwlse variation of the function k^(s ) for 

the delta wing and, inasmuch as the proper osci llat ory lift coefficients 
(for an oscillating airstream.) were not available, the function kgCs) 

was not obtained either. 


SYMBOLS 


A 

b/2 
b(x)/2 

Cl 

c 
c 

c l 
c r 
F(k) 

F 2 (k) 

h(t) 

k = CDCy. J 2 N 


aspect ratio 
semispan of wing 

spanwise coordinate of the leading edge of wing, measured from 
root chord of wing 

steady-state lift coefficient 

section chord 

mean geometric chord 

section lift coefficient 

root chord of wing 

real part of the total lift coefficient for a three-dimensional 
wing oscillating harmonically in pure translational motion 
in incompressible flow 

real part of the local lift coefficient for a three-dimensional 
wing oscillating harmonically in pure translational motion in 
incompressible flow 

vertical velocity of wing 


k-^(s) indicial lift functions for a wing experiencing a sudden change 
in sinking speed, normalized to unity by its steady-state 
value 

kgCs) indicial lift function for a wing penetrating a sharp-edge 
gust, normalized to unity by its steady-state value 

l(x) 

<1 - | PV 2 


lift per unit length 
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k 

S wing area 

s nondimensional distance traveled, root semichords 

t time 

V forward velocity 

Wq amplitude of vertical velocity 

x coordinate along root chord of wing, measured rearward from 

the leading edge 

y spanwise coordinate of wing 

p air density 

to circular frequency 


METHOD 

The method used to obtain the k^_ function was similar to that used 

to obtain the functions for a wing in two-dimensional compressible flow, 
as reported in references 6 and 7. This method, which avoids some of 
the difficulties associated with the direct calculation of the k-j_ func- 
tion, makes use of Garrick's reciprocal relation (ref. 9) between the 
k^ function and the oscillatory lift coefficients for a wing oscillating 

in pure translational motion. 

A survey of sane of the available methods for calculating oscillatory 
lift coefficients indicated that, as a result of limits in ranges of 
aspect ratio and wing geometry implied by various simplifying assumptions, 
no one method should be applied to a.i l the wings considered herein. Con- 
sequently, two methods which appeared to provide a maximum of information 
for a reasonable effort were selected. These two methods, n am e ly, that 
of Lawrence and Gerber (ref. 10) and that of Reissner and Stevens 
(ref. 11), do not directly apply, however, to all wings considered herein 
and, therefore, the consequences of extending the methods beyond their 
original limits were investigated. The results of this study together 
with a discussion of the characteristic features and limitations of each 
method are presented in the sections that follow. 


Method of Reissner and Stevens 

The method of Reissner and Stevens (ref. 11) for the calculation of 
oscillatory lift coefficients is characterized by simplifications which 
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involve a basically two-dimensional analysis of strips taken along the 
span of the wing. The simplifications utilized by Reissner and Stevens 
reduce the surface integral equation to a single -variable integral equa- 
tion in terms of the spanwise variable. The resulting equation is then 
solved by satisfying the boundary conditions along the span of the wing; 
the solution gives values of Fj( k), the real part of the local complex 

oscillatory lift coefficient. The simplifications used in the Reissner 
and Stevens method imply that the method is best suited for wings of high 
aspect ratio; in fact, the limiting two-dimensional case is conta ine d 
directly. 

The accuracy of the results obtained through the use of t he Reissner 
and Stevens method for wings having aspect ratios as low as 3 was inves- 
tigated for elliptical, rectangular, and delta plan forms. Est imat es of 
the accuracy for elliptical and rectangular wings with aspect ratios of 3 
were obtained by comparing k-j_ functions based on the Reissner n-nfl Stevens 

method with a k-j_ function obtained by Robert T. Jones in reference 4 for 
the elliptical wing and with a k^ function based on the method of Lawrence 

and Gerber (ref. 10) for the rectangular wing. The 'results for the ellip- 
tical and rectangular wings are presented in figures 1 and 2 , respectively. 
In both figures the agreement between the curves is considered satisfactory 
for values of s greater than 2. For s less than 2 the results obtained 
by use of the method of Reissner and Stevens and by that of Lawrence »nd 
Gerber are shown dashed to indicate a region where the functions were 
believed to be unreliable for reasons to be discussed in a section to fol- 
low. The accuracy of the results obtained by the Reissner and Stevens 
method for a delta wing having an aspect ratio of 4 was estimated by com- 
paring the Fj(k) function obtained thereby with that obtained from the 

Lawrence and Gerber' method. Appreciable differences in both shape anrl 
magnitude were observed and, inasmuch as the latter method was considered 
to be the more accurate, the Reissner and Stevens method was not used to 
determine either the local or the complete k-^s) function for the delta 

wing. 


Method of Lawrence and Gerber 

The method of Lawrence and Gerber (ref. 10) for the calculation of 
oscillatory lift coefficients is characterized by the reduction of the 
surface integral equation to a single -variable integral equation in te rms 
of a chordwise variable by the choice of a spanwise weighting function. 
The lifting equation is then solved by satisfying the boundary conditions 
at various stations along the chord. The method carries the restrictions 
that the plan form have a straight trailing edge normal to the stream - 
The simplifying assumptions in the method of Lawrence and Gerber imply 
that the method is best suited for wings of low aspect ratio. 
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The accuracy of results obtained through use of the Lawrence and 
Gerber method for wings having aspect ratios as high as 4 was estimated 
by comparing the k-j_ function for a rectangular wing with an aspect ratio 

of 6 based on coefficients calculated by the Lawrence and Gerber method 
with the kj. function obtained by W. Prichard Jones (ref. 5) for the same 

wing. As is indicated in figure 3> in which the two functions are 

presented, there is only a difference of approximately 3 percent between 
the two functions for values of s greater than 2. As in figure 2, the 
function based on coefficients obtained by the Lawrence and Gerber method 
is shown dashed in a region where the function was believed to be 
unreliable . 


Application of Methods 

In view of the results presented in figures 1 to 3* the function 

for the delta wing was calculated from oscillatory lift coefficients 
obtained by use of the method of Lawrence and Gerber; consequently, the 
spanwise variation of local k-j_ functions was not determined for the delta 

wing. Local k^_ functions at various spanwise stations were obtained by 

use of the Reissner and Stevens method, however, for rectangular and, 
elliptical wings. 

Lack of information precluded the calculation of the indicial lift 
function k 2 (s) for gust penetration of wings of finite aspect ratio. 

Although it is shown in reference 9 that the function k 2 (s) can be 
obtained directly from the k-^ function for wings in two-dimensional flow, 

the relation used is one derived from closed-form expressions for the 
lift response to sinusoidal sinking oscillations and to sinusoidal gusts; 
this relation is not valid for wings of finite aspect ratio. The deter- 
mination of the k 2 function requires a knowledge of the lift due to sinu- 
soidal gusts or due to harmonically oscillating flaps for all values of 
flap-chord ratio; this information is not available for wings of finite 
span. The k 2 and k-j_ functions for several wings of vanishingly small 

aspect ratio were determined, however, to aid in the illustration of 
trends of these functions with aspect ratio. A derivation of expressions 
for these functions for vanishingly small aspect ratio is given in the 
appendix, and the indicial lift functions for rectangular-, elliptical-, 
and delta-wing plan forms are presented in table I. 


Method of Calculation 

The indicial lift function for a wing undergoing a sudden change in 
sinking speed can be calculated from the real part of the couple x lift 
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coefficients for a wing oscillating in pure translational, motion "by means 
of the reciprocal relation (ref. 9) 


, / \ 2 r°° F(k) sin ks . . 

k i (s) - “ (1> 


The indicial lift functions k^s) for delta and rectangular "wings were 

calculated from equation (l) "by using the F(k) coefficients presented in 
reference 10 together with additional coefficients which were calculated 
in increments of k = 0.2 by the method of reference 10 for reduced fre- 
quencies k < 2. The validity of the method of reference 10 for very 
high reduced frequencies has been investigated and the method has been 
found to be inapplicable directly for high frequencies. Consequently, 
the F(k) coefficients were faired for 2 < k < « to a value of F(~) = 1. 
This value has been chosen because the method of reference 10 is primarily 
a low-aspect-ratio theory, and for vanishingly Rmal i aspect ratios the 
starting lift k-^(O), which has been shown in reference 7 to be equal 

to F(°°), is presented in the appendix and is equal to 1 . 0 . The fairing 
of F(k) was in the form of a part of a hyperbola, as used in a similar 
problem reported in reference 6. Reference 6 provides a means of esti- 
mating the error due to fairing as it affects the k x function, and for 

the present calculations it is indicated that the error is significant 
only for wing travel s (semichords) less than 2; hence, in this region 
the curves presented are shown by dashed lines . Equation (l) was evalu- 
ated by numerical integration for k < 2 and by analytical integration 
for k > 2 . The kj. functions for rectangular wings with aspect ratios 

of 3 and 6 (the former being presented in fig. 2) are in agreement with 
those derived from reference 11, and it is concluded therefore that the 
F(k) coefficients as calculated from reference 10 are accurate for values 
of k at least as large as 2.0. 

The procedure followed in the calculation of the "local" kj_ functions 
for elliptical and rectangular wings was the same as that for k-j_(s) for 

the delta wing except that the oscillatory coefficients used in equa- 
tion (l) were local values fj(k) as calculated by the method of refer- 
ence 11 , and the Fj(k) curves were faired to the two-dimensional values 
of F-£ (°°) = jt/c L , which was consistent with the method of reference 11. 

The kj_ functions for the complete elliptical and rectangular wings 
were obtained by spanwise graphical integration of the local kj_ functions 



8 


NACA TN 3639 


at six equally spaced stations along the wing semispan. Six stations 
were used rather than the four used in reference 11 because the lift 
distributions obtained from, four stations were believed to be insuf- 
ficiently defined for the purpose of this report. One of the six sta- 
tions chosen was at the wing tip, where the lift was assumed to be zero. 


RESULTS AND DISCUSSION 
Delta-Plan-Form k]_ Function 


The indicial lift functions for a sudden change in sinking Bpeed 
k-|_(s) were determined for delta wings of aspect ratios 0, 2, and 4 and 

are presented in figure 4. It might be noted that, for the delta wing 
as in the cases of wings of other plan forms, a decrease in aspect ratio 
results in a more rapid growth of lift to the steady- state value. 


Indicial Spanwise Loading Distributions 

The local indicial lift functions for a sudden change in sinking 
speed were calculated for six stations along the wing semiBpan on both 
elliptical and rectangular wings of aspect ratios 3 and 6. These data 
for rectangular wings are presented in figure 5 in the form of spanwise 
local lift coefficients per unit wing lift coefficient for several values 
of wing distance traveled. As indicated by figure 5, "the spanwise dis- 
tributions are substantially independent of time, with the possible excep 
tion of the very beginning of motion. No indication of the initial span- 
wise distribution can be given, however, since the theory is not valid 
in this region. The data for the elliptical wings exhibited character- 
istics similar to those for the rectangular wings in that the spanwise 
distributions were, for all practical purposes, equal to the steady-state 
lift distribution, which in this case is elliptical. 


Effect of Plan-Form Shape on the k-]_ Function 

The effects of plan-form shape on the k^ functions for elliptical, 

rectangular, moderately tapered, and delta wings are shown in figure 6. 
The k^ functions for the first three wings of aspect ratio 6 (aspect 

ratio of ^.8k for tapered wing) are grouped together in figure 6(a) 
wherein it is indicated that the k^ functions are almost identical. The 

results for elliptical and rectangular wings of aspect ratio 3, although 
not plotted, show similar characteristics. 
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The k-j_ functions for rectangular and delta wings of aspect ratio k 

are plotted in figure 6(b). Inasmuch, as it has been established that 
the function k^(s) for a rectangular wing is substanti ally the same 

as those for e lli ptical and moderately tapered wings of the same aspect 
ratio, it may be concluded from figure 6(b) that the delta wing exhibits 
a more rapid growth of lift to the steady-state value than do the other 
three plan forms considered. 

The functions k-j_(s) presented in figure 6 were obtained from the 

following sources: The indicial functions for the elliptical wing were 

reproduced from reference k, whereas those for the rectangular and tapered 
wings were reproduced from reference 5 • The delta-wing indicial function 
was reproduced from figure k. 


Effect of Aspect Ratio on the k^ Function 

The effect of aspect ratio on the k-j_ function for delta wings has 

already been indicated in figure 4 and for elliptical, rectangular, and 
moderately tapered wings this effect is indicated in figure 7 , in which 
k-L functions are plotted for rectangular wings having aspect ratios of 0, 

3, k, 6, and 00 . The curves in figure 7 were obtained from the following 
sources : The function for zero aspect ratio was determined by the method 
presented in the appendix of the present report. The functions for aspect 
ratios of 3 and 6 were calculated by means of the method of reference 10. 
The f inactions for aspect ratios of k and °° were obtained from refer- 
ences 5 and 1, respectively. As shown previously, the results for the 
rectangular wings are representative of those for elliptical and moder- 
ately tapered wings. The general effect of decrease In aspect ratio 
for the plan forms considered is a more rapid growth of lift toward the 
steady-state value, as indicated by figures k and 7* 

Inasmuch as the Indicial functions presented herein are normalized 
to unity at the steady-state condition (s = » or k = 0), the applica- 
tion of these functions requires a knowledge of the steady-state lift- 
curve slope. Although the normalized indicial lift functions kj_(s ) 

have been shown to be substanti ally Independent of method of calculation 
(at least for e llip tical, rectangular, and moderately tapered wings), 
the associated steady-state lift-curve slopes can be shown to vary appre- 
ciably with method of calculation. Because a choice of a range of values 
of lift-curve slope for a given normalized function must be made, the 
choice might as well be from a source one considers most reliable. 
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Effect of Aspect Ratio on the kg Function 

The effect of aspect ratio and shape of plan form on the kg function 
is indicated in figure 8 . In figure 8 (a) is plotted the kg function for 

rectangular ■wings of aspect ratios 0 , 4,-and 6 , and in figure 8 (b) is 

plotted the kg function for elliptical wings of aspect ratios 0 , 3 > and 6 . 

The zero -aspect -ratio functions for both wing shapes are those calculated 
by means of the method given in the appendix. The remaining kg functions 

for the rectangular and elliptical wings were reproduced from references 5 
and k, respectively. The effect of aspect ratio is generally the same as 
it was for the k-j_ functions. Although not shown here, it might be inter- 
esting to note that the kg functions for the elliptical and rectangular 

wings of aspect ratio 6 are almost the same, as is the case for the 

k-j_ functions. However, for aspect ratios of 3 or less the kg functions 

differ significantly, so that for low aspect ratios the effect of plan- 
form shape appears to be more significant for the kg function than it is 

for the k|_ function. 


CONCLUSIONS 


The indicia! lift functions for a wing undergoing a change in sinking 
speed in an incompressible medium have been calculated from the oscillatory 
lift coefficients, for a wing oscillating harmonically in pure translational 
motion. These functions are presented for delta wings having aspect ratios 
of 0 , 2 , and K, and the local as well as the total functions are presented 
for elliptical and rectangular wings having aspect ratios of 0 , 3 ; and 6 . 

The present investigation indicates the following results for the 
cases considered herein: 

1. The growth of lift for delta wings is much more rapid than that 
for elliptical and rectangular wings of the same aspect ratio. 

2. For elliptical or rectangular wings, the spanwise distributions 
of the lift are substantially independent of time. 

3 . The normalized unsteady-lift functions are relatively independent 
of plan-form shape for elliptical, rectangular, or moderately tapered 
wings . 
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4. Hie growth, of lift, for all wings investigated, is more rapid 
as the aspect ratio is decreased. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 11, 1956. 
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APPENDIX 

INDICIAL LIFT OF WINGS OF ZERO ASPECT RATIO 


The lift per unit length along a slender body or a wing of very 
low aspect ratio undergoing pure translational motion is shown in refer- 
ence 13 to be 



For pure sinking motion 

h(t) = W Q i(t) (A 2 ) 

where _Z(t ) is the unit jump function. The total lift coefficient can 
then be written as 



Z(x)dx 


(A3) 


where c r is root chord, x* = 2xjcr> 6 = 2Vt/c r , and 6(s) is the 
unit impulse function. 


For a wing penetrating a sharp-edge gust 
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so that 


k 2 (s) 



k2(s) = 1 


where 


*0 


is defined by 



= b. 



These expressions (eqs. (A3) and (A5)) may be evaluated for any 
specific plan form. For instance , the results for rectangular, ellip- 
tical., and delta wings are given in table I. 
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TABLE I 


INDICIAL LI FT FUNCTIONS FOR RECTANGULAR, ELLIPTICAL, AND 
DELTA WINGS OF VANISHINGLY SMALL ASPECT RATIO 


Plan form 

[¥f - [^T 

[ 

^(s) 

k 2 (s) 

Rectangular* 

1 

2( b) + 26(b) 

2(s) 

Elliptical 

1 - (s - l) 2 

i(s) + | 5(s) 

s(2 - b) 0 £ s ^ 1 
2(s) s > 1 

Delta 

s 2 /4 

2(b) + | 6(b) 
5 

s 2 /V 0^6 = 2 

2(s) s > 2 


Vll C\J 
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Method 

“Ref. 1* (Robert T, Jones) 

— Ref. 11 (Reissner and Stevens 
high-aspect-ratio 


Distance traveled, s, root semichords 


Indicial lift functions k-^s) for an elliptical wing of aspect ratio 5 as calculated 
by method of reference 11 and as reproduced from reference 4. 
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Distance traveled, s, root semichords 


Figure 3*- Judicial lift functions kj^s) for a rectangular wing of aspect ratio 6 as calculated 
by method of reference 10 and as reproduced from reference 5. 
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(a) Rectangular plan formj A = 3 . 


(b) Rectangular plan foraj A => 6 . 


ro 

i - 1 


Figure 5 *~ Instantaneous spanwise loading coefficients for several values of distance traveled s. 
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